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ABSTRACT 

Context. The Magellanic Clouds provide a nearby laboratory for metal-poor dwarf galaxies. The low dust abundance enhances 
the penetration of UV photons into the interstellar medium (ISM), resulting in a relatively larger filling factor of the ionized gas. 
Furthermore, there is very likely a molecular gas reservoir probed by the [C n] 157 jum line not traced by CO(l-O), the so-called "dark" 
gas. The Herschel Space Telescope allows us to observe far-infrared (FIR) cooling lines and to examine the physical conditions in the 
gas phases of a low-metallicity environment to unprecedented, small spatial scales. 

Alms. Our objective is to interpret the origin of the diffiise emission of FIR cooling lines in the H n region N IIB in the Large 
Magellanic Cloud. We first investigate the filling factor of the ionized gas. We then constrain the origin of the [C ii] line by comparing 
to tracers of the low-excitation ionized gas and of photodissociation regions (PDRs). 

Methods. We present Herschel/FACS maps of N 1 IB in several tracers, [C n] 157 fim, [O i] 63 fim and 145 jim, [N ii] 122 jim, [N m] 
51 jum, and [O iii] SSjum. Optical images in Ha and [O iii] 5007A were used as complementary data to investigate the eff'ect of 
dust extinction. Observations were interpreted with photoionization models to infer the gas conditions and estimate the ionized gas 
contribution to the [C n] emission. PDRs were probed through polycyclic aromatic hydrocarbons (PAHs) observed with the Spitzer 
Space Telescope. 

Results. [O III] 88 jim is dominated by extended emission from the high-excitation difi\ise ionized gas. This is the brightest FIR line 
throughout N IIB, ~ 4 times brighter than [C ii]. We find that about half of the emission from the ionized gas is extinguished by 
dust. We modeled [O ni] around each O-type star and find that the density of the ISM is < 16 cm~^ on large scales. The extent of the 
[O m] emission suggests that the medium is rather fragmented, allowing far-UV photons to permeates the ISM to scales of > 30 pc. 
Furthermore, by comparing [C ii] with [N ii] 122/im, we find that 95% of [C ii] arises in PDRs, except toward the stellar cluster for 
which as much as 15% could arise in the ionized gas. We find a remarkable correlation between [C ii]+[0 1] and PAH emission, with 
[C ii] dominating the cooling in diff'use PDRs and [O i] dominating in the densest PDRs. The combination of [C ii] and [O i] provides 
a proxy for the total gas cooling. Our results suggest that PAH emission describes better the gas heating in PDRs as compared to the 
total infrared emission. 

Key words. HII regions - photon-dominated region (PDR) - Magellanic Clouds - Galaxies: star formation - Infrared: ISM 



1. Introduction 

Dwarf galaxies are important testbeds for understanding star 
formation in metal-poor environments that are intermediate be- 
tween primordial galaxies and evolved systems such as the 
]V[ilky Way. The low dust abundance in a metal-poor interstel- 
lar medium (ISM) allows far-UV (FUV) photons to permeate 
on larger physical scales in interstellar clouds. This eflfect re- 
duces the size of the CO cores and pushes the C^/C^/CO inter- 
face deeper into the cloud, hence the very low CO luminosities 
(e.g., Leroy et al. 2006) and the difficulty of accurately deter- 
mining the molecular hydrogen reservoir. Since molecular hy- 
drogen (H2) is more efficiently self-shielded from photodissoci- 



Herschel is an ESA space observatory with science instruments 
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ation than is CO, this potentially results in a thick layer where H2 
exists along with C^ and C^ (Grenier et al. 2005; Roellig et al. 
2006; Wolfire et al. 2010). This so-called "dark gas" might rep- 
resent a significant fraction of the molecular mass in metal-poor 
galaxies (IMadden et al. 1997). 

The ratio of the far-IR (FIR) line [C 11] 157 jim to CO(l-O) is 
found to be much higher in dwarf galaxies than in normal galax- 
ies (IVLadden et al. 1997; Pineda et al. 2008), hinting at a possibly 
significant dark gas reservoir. However, many uncertainties re- 
main regarding the exact origin of the [C 11] line as a function of 
star-formation activity and metallicity. Given the ionization po- 
tential of C^ (1 1.3 eV) and the low excitation temperature of the 
157 fim line (91 K), [C n] can originate in the cold and warm neu- 
tral atomic and molecular medium and in the warm ionized gas 
(see also IMadden et al. 1993; Heiles 1994). These phases coexist 
in a single beam for most extragalactic observations, complicat- 
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ing the diagnostics drawn from [C ii] . It is therefore essential to 
probe the smallest physical scale possible in the closest galaxies 
in order to constrain the conditions of [C ii] emission. This has 
proven to be a challenge until now, especially for dwarf galaxies. 

The relatively larger mean free path of FUV photons in low- 
metallicity environments also results in a large filling factor of 
ionized gas. The best probe for estimating this filling factor is 
an optically thin tracer with a low critical density (i.e., the den- 
sity for which the spontaneous radiation rate equals the collision 
rate) and an ionization potential permitting its presence outside 
the dense H ii regions. [O iii] 88yum, [N ii] 122yum, and [N ii] 
205 yum are the best probes available (see critical densities and 
ionization potentials in Table[T]). 

The PACS instrument (Poglitsch et al. 2010) on board the 
Herschel Space Telescope (Pillbratt et al. 2010) observes with a 
spatial resolution of ~ 10-13'', which allows (1) examining the 
distribution of [C ii] in nearby galaxies and (2) comparing the 
local excitation conditions of this important coolant to those of 
complementary tracers of the high-excitation ionized gas ([N iii] 
57//m, [O III] 88yum), of the low-excitation ionized gas (e.g., 
[N ii] 122yum, 205 yum), and of the neutral and molecular gas 
([O i] 63yum, 145 yum). 

The Magellanic Clouds provide both a local dwarf galaxy 
template and a moderately metal-poor environment, with a 
metallicity of 1/2 Z© for the Large Magellanic Cloud (LMC) 
and ^ 1/5 Zo for the Small Magellanic Cloud (SMC). The 
presence of massive star- forming regions, such as SODoradus 
and N 1 1 in the LMC, makes it possible to constrain the ori- 
gin of [C ii] in the various phases involved in star formation. 
The Herschel SHINING Guaranteed Time Key Program (P.I. E. 
Sturm) observed several H ii regions in N 11 with PACS. The 
first objective is to unveil the particular morphology of the ISM 
in a metal-poor environment through the study of the ionized gas 
distribution and filling factor. The second objective is to under- 
stand and constrain the origin of the [C ii] emission in order to 
account for the molecular gas reservoir that is not seen in CO. 
In the present paper, we investigate the large-scale emission-line 
structures in the N 1 IB H ii region that is the northern and bright- 
est H II region in the N 11 complex (Fig.[T]). In a second paper 
(Lebouteiller et al. in preparation), we examine the other H ii 
regions in the Nil complex, with a focus on the physical condi- 
tions in the photodissociation regions (PDRs). 

In Sect.|2] we describe the environment and morphology of 
N 1 IB. We present the observations in Sect.|3]and the data anal- 
ysis in Sect.|4] The spatial distribution of the FIR tracers is ex- 
amined in Sect.|5j and the physical conditions of the ionized gas 
are investigated in Sect.[6] The ionized gas spatial distribution is 
discussed in Sect.lT] Finally, the origin of the [C ii] emission is 
discussed in Sect.[8| We present the Herschel/FACS data reduc- 
tion and analysis tool PACSman in the Appendix. 



2. Environment 

2.1. Morphology 

Nil (also DEM 34; Henize 1956; Davies, Elliot & Meabum 
1976) is the second largest giant H ii region in the Large 
Magellanic Cloud after SODoradus (Kennicutt & Hodge 1986). 
The morphology of N 1 1 consists of nine distinct nebulae 
(NllA-NllI; Rosado et al. 1996) distributed around the OB as- 
sociation LH 9 (NGC 1760), which is dominated by the compact 
cluster HD 32228 (also RadcHff'e 64, Sk-66-28, or Breysacher 9). 
The LH 9 association carved out a cavity and triggered a sec- 
ondary, peripheral starburst (Lucke & Hodge 1970; Parker et al. 



1992; Walborn & Parker 1992). The starburst shell created by 
LH9 reached ~ 120pc in diameter (80x60 pc, Lucke & Hodge 
1970). Figure[T] shows the H i hole corresponding to the cavity 
around LH 9. The nebulae are distributed around the H i shell; 
N 1 IB is the brightest nebula to the north. 

There is evidence of ongoing star formation in N 1 IB. About 
20 HAeBe stars, which are intermediate-mass (3 to 7 M©) pre- 
main sequence stars (1 to 3 Myr), have been detected (Barba et 
al. 2003; Hatano et al. 2006). A methanol maser of 0.3 Jy was 
observed by Ellingsen et al. (1994). Finally, Hatano et al. (2006) 
identified two ultracompact H ii regions, one at (04h56m48s, 
-66° 24' 34''), corresponding to the radio source known as W2 
(or B0456-6629) in Indebetouw et al. (2004), and the other at 
(04h56m57s -66°25'13"), which is known as Wl (or B0456- 
6629). 

From the Hubble Space Telescope (HST) Advanced Camera 
for Surveys (ACS) archival images (Fig. [2]), two regions can 
be distinguished within N 1 IB with distinct morphological at- 
tributes. The eastern half has sharp edges toward the east and the 
south, indicating the presence of relatively dense interfaces. A 
lower density zone is seen in the middle of N 1 IB, with the OB 
association PGMW 3204/09 (hereafter the "footprint" cluster; 
Sect. 2.2) surrounded by several small clouds with bright ion- 
ization fronts, including the kiwi-shaped cloud and the Y-shaped 
cloud described in Naze et al. (2001). The western half of N 1 IB 
is dominated by PDR structures north of the local OB associa- 
tion LH 10. Several ionization fronts (II, 12, and 13 in Barba et 
al. 2003) are seen above LH 10. 

In the following we concentrate on four distinct regions that 
were selected based on their extreme nature among N 1 IB, trans- 
lating into significant diff'erences in line ratio diagnostics: 

- The fragmented region within the orange box in Fig.|2j which 
we refer to as W2 from now on, which contains relatively dense 
interfaces under the influence of nearby massive stars from the 
LH 10 cluster. 

- The stellar cluster below W2, which we refer to from now on 
as LH 10 (red box). LH 10 is embedded in a relatively difl'use 
gas. 

- The northern region (green box), where the absence of optical 
edge or arcs suggests a smooth transition into the general, lower 
density, LMC environment. This region is little, if at all, afl'ected 
by known massive stars. 

- The compact region Wl (blue box), which contains a small 
stellar cluster surrounded by dense shells. 



2.2. Main ionizing sources 

The stellar cluster associated with N IIB is the LH 10 OB asso- 
ciation (also NGC 1763, IC 2115, Lucke & Hodge 1970). LH 10 
is the youngest cluster in N 11 (< 3 Myr; Walborn et al. 1999). 
Most of the stars in LH 10 are gathered into several compact 
complexes located south of the W2 star-forming region. The 
PGMW 3128 and the PGMW 3120 complexes are located near 
W2, while the PGMW 3090 and PGMW 3070 are farther away 
to the south (Parker et al. 1992; Fig. [2]). From the high-resolution 
HST/ACS image, we count at least 260 stars south of W2. 
Including faint sources increases the number to more than 300 
stars. This is a lower limit since some complexes are unresolved. 
At least two 03 stars were identified in the southern LH 10 clus- 
ter, PGMW 3058 and PGMW 3061, sHghtly off' the area covered 
by our PACS observations (Parker et al. 1992). The presence of 
"intact" (not yet evolved) 03 stars suggests that no supernova 
event has occurred yet. 
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Fig. 1. The giant H ii region Nil, with Spitzer/MIFS 24/jm shown in red, Spitzer/IRAC S.Oyum in green, and H i column density 
(from Kim et al. 2003) in blue. The Spitzer/MIFS 24 jim emission is dominated by very small grains that are stochastically heated 
and by big, warm, dust grains in thermal equilibrium with the interstellar radiation field. Spitzer/IRAC 8.0 yum emission is dominated 
by PAHs. The polygon shows the Herschel/FACS spectroscopic observations of N IIB presented in this paper. 



The "footprint" cluster is located to the east of NllB 
(see lower left panel in Fig. [2]) and close to the kiwi-shaped 
cloud, with several massive stars and concentrations identified 
by Parker et al. (1992). The main sources are PGMW3204- 
3209-321 1, with at least 30 stars (Fig.[2]). An 03 III star was dis- 
covered in PGMW3209 (star 'A" in Walbom & Parker 1992). 
Another concentration is located nearby, to the east of the Y- 
shaped cloud with at least a dozen stars. 

2.3. Electron density 

Mac Low et al. (1998) find an electron density of ^ 38 cm"^ 
from the Ha emission measure, assuming N IIB is a uniform 
ionized sphere. By estimating the physical size of structures 
within the nebula. Naze et al. (2001) find that the emission mea- 
sure is consistent with ~ 15 cm"^ in the lowest surface bright- 
ness regions. It reaches 50 - 150 cm"^ in moderately bright arcs 
and up to 200-500 cm"^ in the bright ridges of dust clouds (e.g., 
the kiwi- shape cloud). The bright optical arcs studied by Naze et 
al. (2001) correspond to the highest densities and also show low 



[O III] 5007A/HQf ratios, which is attributed to lower ionization 
and excitation. 

Tsamis et al. (2001) observed the [S ii] A6131/A6116 and 
[O ii] A3129 / A3126 doublets in long-slit spectroscopy (see slit 
position in Fig.|2]) and find rig = 80 cm"^ and llOcm"^, re- 
spectively. They measure a significantly higher density from the 
[CI III] .15537//15517 doublet, with 1700 cm'^ which they at- 
tribute to the higher critical densities involved in the [CI iii] 
transitions. Strong density variations are therefore expected in 
N 1 IB. We tentatively expect a density of ~ 25 cm"^ for the low- 
density nebula and densities between 100 - 1000 cm"^ for the 
bright optical arcs. 

2.4. Chemical abundances 

Chemical abundances in NllB were derived by Tsamis et al. 
(2003), who find 12 -h log(0/H) = 8.41 and 12 -h log(N/H) = 
6.93. Using the oxygen solar abundance 12 -h log(0/H) = 8.69 
from Asplund (2009), the metalHcity is 0.527©. The N/0 abun- 
dance ratio is significantly lower than the solar value since N/H 



3 



Lebouteiller et al.: Physical conditions in the gas phases of the H ii region LMC-N 1 IB unveiled by Herschel 




15.0 10.0 05.0 4:57:00.0 55.0 50.0 45.0 40.0 35.0 56:30.0 

Right ascension 




Fig. 2. Top - HST/ACS color-composite image of the NllB region (blue: F814; green: F658N/Ha; red: F660N/[N ii]). 
Images were downloaded from the Hubble Legacy Archive {http://hla.stsci.edu/) and later stitched with Montage 
(http://montage.ipac.caltech.edu/). The white stripes correspond to a gap between individual exposures. The yellow polygon in- 
dicates the field-of-view of the Herschel/FACS observations 2.3' x 2.3'). The 4 rectangles in color indicate specific regions 
discussed in the text. The dashed circles represent the ISO-LWS observation (from 65" to 80" FWHM beam) and the slit represents 
the optical observation by Tsamis et al. (2001). Bottom - Close-up of the eastern (left panel) and western (right panel) sides. The 
stellar complexes are labeled as the brightest star with a '-h' sign appended. 
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is 7 times below the solar abundance. This nitrogen deficiency 
is observed in several H n regions of the LMC and is likely due 
to depletion of N on dust grains (e.g., Garnett 1999). For the car- 
bon abundance, we assume 12 + log(C/H) = 7.9 which is the 
average abundance in LMC H ii regions (Garnett 1999). 

3. Observations 

N IIB was observed on January 4, 2010 with the PACS spec- 
trometer on board Herschel. The main FIR fine-structure lines 
were mapped in the wavelength- switching mode (later decom- 
missioned and replaced by the unchopped scan mode): [N iii] 
57yum, [O i] 63yum, [N ii] 122 yum, [O i] 145 yum, [C ii] 157 yum, 
[N II] 205 yum (OBSID 1342188940), and [O iii] 88yum (OBSID 
1342188941). Two follow-up observations were performed in 
unchopped scan mode, one pointing toward W2 in [N ii] 122 yum 
on April 4, 201 1 (OBSID 1342219439), and one pointing toward 
Wl in [N II] 122 yum and [O i] 145 yum on July 19, 201 1 (OBSID 
1342225175). The [O iii] 52yum was not observed with PACS 
because of the low transmission at this wavelength, but the line 
was detected with ISO (Sect. |4.2| ). The observed lines are listed 
in Table[T] along with a summary of the main transition parame- 
ters. From now on we use the line notation with the wavelength 
as subscript for ambiguous cases (e.g., [O i]63 for [O i] 63 yum). 

As explained in Poglitsch et al. (2010), the PACS spectrom- 
eter is an integral-field spectrometer that samples the spatial di- 
rection with 25 pixels and the spectral direction with 16 pixels. 
Each spectral pixel "sees" a distinct wavelength range that is 
scanned by varying the grating angle. The resulting projection 
of the PACS array on the sky is a footprint of 5 x 5 spatial pix- 
els ("spaxels"), corresponding to a 47'' x 47'' field-of-view. At 
a distance of the LMC (48.5 kpc), the footprint size 47" corre- 
sponds to 1 1 pc. The map observations of N 1 IB were designed 
as 3 X 3 raster map, with an overlap of 1.5 spaxels. Figure|2] 
shows the PACS map coverage (about 2.3' x 2.3') on the opti- 
cal image. According to the PACS Observers ManuaQ the point 
spread function full width at half maximum (FWHM) is ^ 9.5" 

2.2 pc) between 55 yum and 100 yum, and it increases to about 
^ 13" (^3pc)atl80yum. 

4. Data analysis 

4.1. Herschel/PACS 

The data reduction was performed in HIPE track number 
4.0.700 (Ott 2010) using a combination of the default 
wavelength- switching and unchopped scan scripts. 

4.1.1. Wavelength-switching 

A spatially off'set position is usually observed at the beginning 
and at the end of the observation to subtract the telescope back- 
ground emission. In our mapping observations, the off'set was 
located on the source itself (central raster position in the map) 
so that it was not subtracted but instead used to increase the 
signal-to-noise ratio (S/N) by ~ V2 at that location. Based on 
the uncho pped sc an observation, for which an off'set was avail- 
able (Sect. |4.1.2| ), we estimate that the lack of an off'set in the 
wavelength- switching observation resulted in a lower S/N in the 
final spectra. No line contamination from the Milky Way was ob- 
served in the off'set in the unchopped scan observation, implying 
that our line measurements are not aff'ected by foreground line 

^ http .'//herschel. esac. esa. int/Docs/PACS/html/pacs-om. html 



emission, as expected from the low galactic latitude of LMC- 
Nll (~ -66.4°). 

The wavelength- switching mode was initially designed to re- 
move the variations in the telescope background emission on 
small timescales by switching wavelengths rapidly, resulting in 
a diff'erential spectral line profile. The wavelength modulation of 
the grating follows the AABBBBAA pattern, where A observes 
at the nominal wavelength and B observes the modulated wave- 
length (with a shift of about half of the FWHM). 

The diff'erential line profile resulting from the wavelength- 
switching proved itself difficult to fit, especially for low S/Ns. 
We decided to consider the line profile in the combined spectrum 
A-hB (i.e., without subtracting the switched wavelength spec- 
trum). This is eff'ectively a similar approach to the unchopped 
scan mode. The final spectrum was divided by a factor 1 .3 for the 
blue camera and 1.1 for the red camera in order to account for 
diff'erences between ground and in-flight performances (PACS 
ICC calibration document PICC-KL-TN-041). 



4.1.2. Unchopped scan 

Follow-up pointed observations of [N ii]i22 and [O i]i45 were 
performed with the unchopped scan mode. The ofl'set exposures 
were checked for line emission and were averaged before sub- 
traction from the "on- source" observations. The unchopped scan 
observations took longer (per line and for one raster position) 
than the wavelength-switching, and transient cosmic rays were 
observed on the timeline of several photoconductor pixels. We 
corrected the signal by using a multiresolution algorithm that is 
part of the data reduction and analysis P AC Smart (see Appendix). 

4.1 .3. Line fitting and map projection 

Line fits were performed by PAC Smart (see Appendix) for each 
spaxel at each raster position. A Gaussian profile was adjusted 
simultaneously with a polynomial continuum of degrees one 
to three. All lines have heliocentric radial velocities of 305 ± 
20 kms"^ with no significant variations across the map. This 
velocity is consistent with the value 292 ± 23 kms"^ found by 
Rosado et al. (1995) for the entire N 1 IB nebula. The instrument 
spectral resolution ranges from ^ 55 kms"^ to ^ 320 km s"^ de- 
pending on the band/order and wavelength. No significant ad- 
ditional broadening is observed for the long-wavelength lines 
[O i]i45 and [C ii], while a broadening of 70 - 90 kms"^ is ob- 
served for [N III], [O i]63, and [O iiilgg (for which the spectral 
resolution is finer), which might be real. Example of line fits are 
shown in Fig. [3] 

For the mapping observations, the line fluxes of all the spax- 
els are projected on a subpixel grid (see Appendix). The subpixel 
size (1/3 of the spaxel size, i.e., ^ 3'' or ^ 0.7 pc at the distance 
of the LMC) is chosen in order to recover the best possible spa- 
tial sampling in the footprint overlap regions. For comparison, 
the map of the Orion bar in Bernard-Salas et al. (201 1) is 0.25 pc 
in size, i.e., a fraction of one pixel in our maps. 

The S/N is > 50 in the [O miss fine across the map, > 12 
in the [C ii] map, between 1 and 100 in the [O i]63 map, be- 
tween and 17 in the [O i]i45 map, and between 0.5 and 7 in 
the [N III] map. The [N ii]i22 fine is barely detected in a few 
positions, while [N ii]205 is not detected anywhere because of 
the known spectral leakage in the Rl band (see PACS Observers 
Manual). Final maps are shown in Fig.[4j 
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Table 1. Main parameters of the transitions used in this study. 



Line 


A 


IPi^ 


IP2^ 


Configuration 




d 








(eV) 


(eV) 




/ -3\ 

(cm ^) 


/ -3\ 

(cm ^) 


(K) 


[0 l]63 


63 jim 




13.62 


3d 3d 




^ (4 - 10) X 10 


228 


IJl45 






i J.DZ 


3 p 3 p 




~ i^i — ZU^ X iU 


J A 1 


[C II] 


157//m 


11.26 


24.38 


^^3/2 - 


^ 80 


^ 3000 - 5000 


91 


[0 III]52 


52 Jim 


35.12 


54.93 


'Pi -'Pi 


^ 3 500 




441 


[0 III] 88 


88yL/m 


35.12 


54.93 


'Pi - 'Po 


^510 




163 


[N III] 


57 //m 


29.60 


47.45 


^Ps/l - ^Pl/2 


^ 3000 




251 


[N II] 122 


\22/jm 


14.53 


29.60 




^400 




188 


[N Il]205 


205 yum 


14.53 


29.60 


^Pi - ^Po 


^ 180 




70 


[0 IIl]opt 


5007A 


35.12 


54.93 




^7x10^ 




2.9 X 10^ 



Ionization potential to create the species. Ionization potential to ionize the species. Critical density for collisions with electrons; i.e., the 
density for which the spontaneous radiation rate equals the collision rate. Critical density for collisions with H i and H2. 



Table 2. Line fluxes derived from the ISO LWS spectrum. 



Line 


ISO^ 




ISO^ 




PACS^ 




[OlIl]52 


1.2x 10" 


13 


8.4 X 10- 


14 






[N III] 


1.5 X 10" 


14 


1.1 X 10- 


14 


9.9 + 5.4 X 10- 


-15 


[Ol]63 


2.7 X 10- 


14 


1.9 X 10- 


14 


2.1+0.1 X 10- 


-14 


[O III] 88 


1.7 X 10- 


13 


8.5 X 10- 


14 


9.8 + 1.7 X 10- 


-14 


[Nil] 122 


< 9.1 X 10-16 


< 2.8 X 10-16 


6.1+5.1 X 10- 


-16 


[0 I] 145 


1.2 X 10- 


15 


4.2 X 10- 


16 


1.1 + 0.3x10- 


-15 


[C II] 


2.9 X 10- 


14 


1.0 X 10- 


14 


2.7 + 0.7 X 10- 


-14 



Notes. Line flux is given in Wm-^. 

(«) Point-like source calibration. Errors are ~ 25% for [O i]i45 and 
~ 10% for the other lines. The effective aperture radius is ^ 35'' for 
[O i]i45 and [C 11], ^ 39'' for [N ii]i22, ^ 40" for [O iii]88, ^ 42" for 
[O iii]52 and [N iii], ^ 43.5" for [O i]63. Extended source calibration. 
Integrated PACS flux within the corresponding ISO aperture. 



4.2. Ancillary data 

N 1 IB was observed with the ISO space telescope (Kessler et 
al. 1996) in the LWS filter (45 - 200 yum) with a beam size be- 
tween 65'' and SC. We retrieved the archival data, extracted the 
spectrum (Fig.[5|, and measured the line fluxes. The flux deter- 
minations are listed in Table[2] The values assuming a point-like 
source calibration agree within 10% with Vermeij et al. (2002^ 
We also extracted the line fluxes assuming an extended source 
calibration. The good alignment of the dust continuum across 
the spectral orders in Fig. [5] suggests that the dust distribution 
is point-like. We compared the ISO line fluxes with the inte- 
grated fluxes derived from the PACS maps in the same aperture 
(Fig.|2]). The PACS fluxes agree remarkably well with the ISO 
results. The PACS fluxes are slightly underestimated because of 
the incomplete coverage of the ISO beam (Fig.|2]). 

We also used the calibrated HST/WFPC2 images in the 
F502N and F658N narrow band filters, kindly provided to us by 
Y. Naze (see also Naze et al. 2001). The F502N band is domi- 
nated by the [O iii] 5007 A (hereafter [O iii]opt) fine emission. The 
continuum emission, which is dominated by dust- scattered light 
in N 1 IB, represents only a small fraction (< 1%) of the emission 
within the narrow filter as compared to [O iii]opt in the long-slit 
used by Tsamis et al. (2008). We assumed that the continuum 
emission remains negligible throughout the region. A similar as- 
sumption was used for the F658N band and the Ha line. The 
WFPC2 images were convolved with a Gaussian PSF to reach a 
spatial resolution of 12'', i.e., slightly larger than the PACS reso- 



lution (^ 10''). The WFPC2 images were then resampled to the 
PACS image pixel size (Sect. 4.1.3] ). The final [O iii]opt and Ha 
maps are shown in Fig. [6] 

We derived the polycyclic aromatic hydrocarbon (PAH) 
emission by using the archival Spitzer/IRAC images. We used 
band 4, centered at S.Oyum, which observes the PAH complex 
at 1.1 /im and 8.6yum along with the warm dust continuum. To 
remove the warm dust contribution, we assumed a power-law 
emission on the Wien side of the spectral energy distribution 
(SED) and used the emission in bands 1 (3.6 jim) and 2 (4.5 yum) 
to extrapolate to the band 4 wavelength range for each pixel. 
The correction thereby applied is ~ 10%. We then calibrated the 
flux density in band 4 by using a template source spectrum, as 
explained in the Spitzer synthetic photometry cookbool|^ The 
template spectrum was built from an incomplete Spitzer/IRS 
map of N IIB, covering the LH 10 and W2 regions (AORKEY 
22469632). From the IRAC transmission curve, we further es- 
timated that about 45% of the total PAH emission (including 
features from 5 yum to 18 jim) falls within band 4. We therefore 
corrected the PAH map by this factor. The final PAH map is pre- 
sented in Fig. [6] 



5. Distribution of the FIR lines 

The main observational results (noise level and integrated flux) 
are provided in Table[3] We discuss the detailed results for each 
line in the following. 



5.1. [Oi] 

Both [O i]63 and [O i]i45 have high critical densities (Table[T]) and 
trace mostly the dense gas in PDRs. The two lines show nearly 
identical distributions (Fig.|4]), most of the difl'erences being due 
to the lower S/N of the [O i]i45 map in the fainter and more 
difl'use regions. The [O i]i45 emission is dominated by the region 
W2, while fainter emission is also seen toward the three compact 
regions on the eastern side, including the region Wl. 

The [O i]63 fine is detected throughout the entire map, its 
intensity varying by a factor of ~ 90 (TablejSj. The line ratio 
[O i]i45/[0 i]63 varies from ^ 0.05 to 0.11 across the map. 
The integrated [O i]63 emission in the map is 3.9 x 10"^^ Wm"^, 
about ^ 39% of which comes from W2, and ^ 44% come from 
both Wl and W2. This fraction is remarkably large as compared 
to the area covered by W1+W2 within the PACS map (~ 13%). 



2 As noted in Vermeij et al. (2002), the "N 1 1 A" LWS observation is 
in fact pointed toward N 1 IB. 



http://irsa. ipac. caltech. edu/data/SPITZER/docs/dataanalysistools/ 
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Fig. 3. Line fits are shown for one spaxel observing the W2 region. The histogram shows the rebinned data, the thin curve shows the 
continuum, and the thick curve shows the fit. The fit was performed on the data cloud, which contains ^ 140000 points. We show 
the rebinned spectrum for display purposes. 



5.2. [Cii] 

The [C ii] line is thought to be the dominant cooling mechanism 
in low-density PDRs and in the diff'use ISM for temperatures 
< 8 X 10^ K (Dalgarno & McCray 1972; Tielens & Hollenbach 
1985a, 1985b; Wolfire et al. 1995). Depending on the ionization 
fraction of the gas, excitation of the ^^3/2 level is due to inelastic 
collisions with neutral hydrogen atoms and molecules or with 
electrons (Dalgarno & McCray 1972; Stacey 1985; Kulkarni & 



Heiles 1987). Higher electronic states (^P, ^D, and ^S) lie more 
than 6.2 x 10"^ K above the ground state, temperatures that are 
not attained in most astrophysical environments (see Crawford 
et al. 1983). 

[C 11] varies by a factor of 20 across the map (Fig.|4]). 
Although it shows a flatter distribution than [O i]63, [C 11] is 
the brightest where [O i]63 is the brightest, i.e., toward W2. 
The emission around W2 is, however, more extended than the 
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Fig. 4. HerschelfPACS maps (in J2000 c oord inates). The colorbar gives the flux in Wm ^ sr ^. The dashed rectangles correspond 
to the zones described in Fig.|2]and Sect. 2.1 The spatial resolution is about ^ 9.5'' for [N iiilsv, [O i]63, and [O iiijss, ~ 11.5'' for 
[N II] 122, [O i]i45, and [C iijis?. 



[O i]63 emission, following the main optical arc II well (Fig. [2]). 
The [O i]63/[C ii] ratio shows strong variations, from 1.5-2 
toward Wl and W2 down to 0.3 everywhere else. The total 
[C ii] emission in the map is 8.4 x 10"^"^ Wm"^, about ^ 19% 
of which comes from W2. This fraction is only slightly greater 
than the area fraction covered by W1+W2 within the PACS map 



(~ 13%), suggesting that most of the [C ii] emission is located 
in relatively extended and difl'use regions. 

The [C ii] line was observed in N 1 IB with the NASA Kuiper 
Airborne Observatory (KAO) by Israel & Maloney (2011), who 
find a peak emission of 2.2 x 10"^Wm"^sr"^ with a beam 
of ~ 68'', and by Boreiko & Betz (1991) who found 1.0 x 
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Fig. 5. ISO LWS spectrum of N 1 IB. See Fig.|2|for the location of the observation pointing. The colors correspond to the different 
spectral orders. The [O miss is observed in 2 spectral orders with different spectral resolutions. The displayed spectrum was extracted 
using a point-like source calibration. 



Table 3. Summary of observations and integrated emission prop- 
erties. 



Line 




RMS^ 
(Wm-2) 




SNR" 


^tot 

(Wm-2) 


[Ol]63 


Map 


2.1 X 10- 


17 


0/6/86 


3.9 X 10-^^ 


[0 I] 145 


Map 


4.7 X 10- 


18 


0/2/35 


> 3.4 X 10-15 




Pointing 


1.5 X 10- 


18 


0/5/45 




[C II] 


Map 


9.7 X 10- 


18 


10/43/95 


8.4 X 10-1^ 


[0 III] 88 


Map 


2.0 X 10- 


17 


29/101/152 


3.8 X 10-13 


[N III] 


Map 


3.0 X 10- 


17 


0.7/4/14 


> 3.2 X 10-14 


[N II] 122 


Map 


8.6 X 10- 


18 


0/1.4/5 


> 3.2 X 10-15 




Pointing 


1.4 X 10- 


18 


0.3/5/13 




[N Il]205 


Map 


5.5 X 10- 


18 


0/0.5/2 


> 1.0 X 10-15 



Radial velocity. Median root mean square (RMS) value per 
spaxel across the map (or across the footprint for the pointed 
observations). Minimum, median, and maximum S/N ratio per 
spaxel. Integrated flux in map. 



IQ-^Wm-^sr"^ in a 43'' beam. With BICE, IMochizuki et al. 
(1994) measured 6 x 10"^ Wm"^ sr"^ in a 12.4' beam. The peak 
emission we detect (8x10"^ Wm"^ sr"^ with a beam of ~ 11.5") 
is significantly larger by a factor of ~ 3.5 than in Israel & 
IVLaloney (2011). We attribute this difference to the relatively 
smaller beam of PACS (^ 12"; Sect.[3]). 



5.3. High-excitation ionized gas tracers 

Both [O III] 88 and [N iii] trace the gas photoionized by massive 
stars. In the N 1 IB region, we identify at least two 03 stars, and 
probably as many as a fe w do zen later-type O stars able to create 
and ions (Sect.[l2|). 

The critical density of [O iii]88 is ^cr ~ 510cm"^ (as com- 
pared to ^ 3000 cm"^ for [N iii]), which makes it the best tracer 
available for low-density, high-excitation, ionized gas. The elec- 
tron density in N 1 IB is expected to be around 25 cm"^ on a 
global scale, and reaches up to a fe w hu ndred per cm^ in low- 
ionization bright optical arcs (Sect. |2.3| ). We therefore expect 
most of the gas to be able to de-excite radiatively through both 
[O III] 88 and [N iii]. 

Both [O III] 88 and [N iii] show a fairly similar spatial distribu- 
tion with a relatively narrow range of values (both within a factor 
of 7; Fig.|4]). IMost of the emission is concentrated between W2 
and the stellar cluster LH 10 to the west, and around the "foot- 
print" stellar cluster to the east. The total [O iii]88 emission in 
the map is 3.8 x 10"^^ Wm'^, ^ 14% of which comes from W2, 
which is only barely larger than the actual area covered by W2 
within the PACS map. [O iii]88 is therefore dominated by ex- 
tended emission, making it the brightest FIR line over the entire 
map area. [O iii]88 is brighter than [C ii] by a factor of 1.2 to 20 
throughout the map, and its integrated flux is four times brighter 
than [C ii]. The predominance of [O iii]88 as compared to [C ii] 
is also found in the pointed ISO observatio n, wi th [O iii]88 being 
about seven times brighter than [C ii] (Sect.|4.2[). 
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Fig. 6. [O iiiJopt, Hof, and PAH maps, from top 
images were convolved at the resolution of the 
map. See Fig. |4] for the figure description. 
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to bottom. The 
Herschel [C ii] 



The peak of [O miss lies directly south of W2, where it is 
slightly off'set with respect to the location of the main ionizing 
front "II" (Fig. [2]). The [N iii] peak lies between the [O iiilgg peak 
and the PDRs. Considering the ionization potentials of and 

(35.1 eV and 29.6 eV respectively), the diff'erence in the peak 
location suggests that the ionization structure on the western half 
of N 1 IB is dominated by the LH 10 cluster in the south. In the 



eastern half, the emission is concentrated toward Wl, the "foot- 
print" cluster, and the optically bright filaments (Fig. [2]). 

In the following (Sect. |6.3| ) we compare [O iiilgg to the opti- 
cal [O iiiJopt fine at 5007A (Fig. [6]). The [O iiiJopt fine has a much 
higher critical density of ^cr ~ 8.6 x 10^ cm"^ 510cm"^ for 
[O miss), and thus traces a wider range of gas densities, pro- 
vided there is no significant extinction by dust. The [O 

III] opt a^nd 

[O III] 88 fines also have diff'erent excitation temperatures. Their 
theoretical ratio is shown in Fig. [7] as a function of density and 
temperature (see also Osterbrock & Ferland 2006 and Palay et 
al. 2012). Tsamis et al. (2001) found ~ 9400 K by observing 
the [O III] optical lines /14363, /14959, and /15007 in a long slit. 
For such high temperatures, the [O 

iii]88/opt ratio mostly depends 

on density. 

The [O III] 52 line was not observed with PACS, but it 
was observed and detected with ISO (Sect. |4.2| ). The observed 
[O III] 88/52 ratio lies between 1 and 1.4. The exact value depends 
on the asumed flux calibration (extended or point-like). Within 
the ISO beam, we estimate that only ~ 20% of the [O iii]88 flux 
comes from W2, while the rest comes from extended emission. 
We thus consider that the [O m] 88/52 ratio should be close to 
the extended source calibration value, i.e., ^ 1. Figure]?] shows 
that such a ratio implies an electron density of around 100 cm"^. 
This value is somewhat higher than the density derived from the 
emission measure in the global ISM of N IIB (Sect. |2.3] ). We 
note, however, that 100 cm"^ is an upper limit since the [O iii]88 
emission is not completely uniform. 



5.4. Low-excitation ionized gas tracers 

The [N ii]i22 and [N ii]205 fines trace the low-excitation ionized 
gas (the ion exists for energies > 14.5 eV). The critical den- 
sity of [N ii]i22 (^cr ~ 400 cm"^) is similar to that of [O iii]88 
(Table[T]). [N ii]i22 was only barely detected toward a few posi- 
tions in the mapping observation because of a too short integra- 
tion time. The upper limit on the total emission in the map is 

< 5.5 X 10"^^ Wm"^. The follow-up pointed observations suc- 
cessfully detected [N ii]i22 toward W2 and Wl, with peak emis- 
sions of 1.4x 10"^ Wm"^ sr"^ and 8.5 x 10"^ Wm"^ sr"\ respec- 
tively. 

Unfortunately, [N ii]205 was not detected in the mapping ob- 
servation. The upper limit on the total emission in the map is 

< 3.6xlO"^^W m"^. Since the [N ii] 122 is systematically brighter 
than [N ii]205 for densities over 10 cm"^ (Oberst et al. 2006), the 
upper limit on [N 11] 205 is not a useful constraint on the gas den- 
sity and on the models discussed in this paper. 



6. Physical conditions of the ionized gas 

6.1. Contribution from shocks 

The presence of shocks from supernovae remnants or from stel- 
lar winds could modify the electronic level population and there- 
fore line ratios. A difl'use soft X-ray emission component was 
detected in N IIB (Mac Low et al. 1998; Naze et al. 2004), and 
wind-blown bubbles around individual stars or groups of stars 
were detected kinematically by Naze et al. (2001). These bub- 
bles are all X-ray emitters despite their low expansion veloci- 
ties of 10 - 15kms"^. On the other hand, as noted in Naze et 
al. (2001), the young age of the LH 10 association, as opposed 
to the central association of N 1 1 (LH 9), suggests that the most 
massive stars did not have time to evolve into supernovae. Strong 
interstellar shocks are therefore not expected in N 1 IB. 
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Table 4. Main parameters for the models. 




0.0 



12 3 4 

log density [cm-3] 

Fig. 7. Density diagnostic plot using the [O iiilss/opt ratio (top) 
and the [O m] 88/52 ratio (bottom). The [O iii]88/opt ratio depends 
mostly on the electron temperature for T < 8000 K. Results were 
computed with Cloudy using a half-solar metallicity (Sect.|6.2). 



To verify these assumptions, we compared the observed line 
ratio [O iii]88/opt to the predictions from the shock model de- 
scribed in Allen et al. (2008). For LMC abundances and a den- 
sity of 1 cm"^, the [O 

iiiJss/opt ratio in the shocked gas is always 
lower than ^ 0.23 throughout the range of velocity and mag- 
netic field values in the model. A somewhat lower ratio would 
be obtained for densities higher than 1 cm"^ . Accounting for the 
extinction by dust (which aff'ects the optical line emission at least 
toward W2 and LH 10; Sect. |6.3| ), our observations show that 
[O iii]88/opt ^ 0-25 across the nebula (Fig.lSb). As a result one 
can see that (1) shocks will tend to lower the observed [O iii]88/opt 
ratio, and (2) although shocks might exist, they do not dominate 
the ionization structure of the N 1 IB nebula. 



6.2. Photoionization modeling 

In the following, we model the physical conditions in the N 1 IB 
nebula with the photoionization code Cloudy (c 10.00; Ferland et 
al. 1998). A single stellar source is used to model the radiation 
field. We used the stellar atmosphere grid WM-Basic (Pauldrach 
et al. 2001 and references therein) because it includes non-LTE 
eff'ects. The stellar grid was computed with half- solar metallic- 



Parameter 


Value 


Gaseous abundances'^ 


0:8.40, N: 6.94, C: 7.91 


Ionization parameter U 


[-5,0] 


Stellar temperature 


[30000, 50000] K 


Stellar atmosphere metallicity 


0.5 Zo 


Density (constant) 


[1,10^] cm-3 


Filling factor 


1 



Notes. Chemical abundances, expressed as 12+log(X/H), are discussed 
in Sect.Ea 



ity and main- sequence parameters. The intensity of the radiation 
field varies through the ionization parameter U, which is defined 
as 



U = 



0(H) 
n(B)c' 



(1) 



where 0(H) is the flux of ionizing photons per unit surface and 
^(H) is the hydrogen density. The other varying parameters are 
the stellar radiation efl'ective temperature Tq^ and the gas density 
n (assumed to be constant). The main parameters for the models 
are listed in Table|4j IModels are stopped at the ionization front, 
i.e., when the electron fraction reaches 0.5. Figure[8] shows the 
line ratios predicted by the models, together with the observa- 
tions toward specific regions in the map (the compact regions 
Wl and W2, the stell ar cl uster LH 10, and the relatively diff'use 
northern region; Sect. |2.1| ). 

The observed [N iii]/[0 iii]88 ratio is compatible with ioniza- 
tion parameters > -3 and with stellar temperatures between ~ 
40000 K and ~ 45000 K (Fig.[8^). Such temperatures correspond 
to 07-05 main sequence stars, many of which are detected in 
N IIB (Sect.jT]). Only a couple of more massive stars than 05 
were detected that must drive T^^ slightly higher. Inversely, radi- 
ation temperatures colder than 40000 K fail to reproduce the ob- 
served [N iii]/[0 III] 88 ratio because corresponding models pro- 
duce too little [O III] 88 emission. 

Considering that the gaseous nebula in N II B should have 
a density higher than > 10 cm"^ (Sect. 23), the average 
[N iii]/[0 III] 88 ratio in each region (Fig.[8^) is compatible with an 

ionization parameter 2 for a radiation temperature of T^^ = 

40000 K and 3 for T^^ = 45000 K. Stifl, [N iii]/[0 iii]88 

shows significant variations across NllB, by a factor of ~ 5. 
The lowest value is found in LH 10 and agrees with both low- 
density and high ionization-parameter. Inversely, the highest 
value is found in the northern region and it agrees with a low 

ionization-parameter ( 3), which is expected from the larger 

distance to the known ionizing clusters. 

For illustration, an ionization parameter of -2 would corre- 
spond to a 25 cm"^ cloud located 6 pc from a 05 V star (U = -1 
implies a distance of 2pc while U = -3 implies 18pc) or 
to a 200 cm"^ cloud located 2pc away. A lower limit of 5pc 
can be set on the distance between the W2 ionization fronts 
and most stars in LH 10 by measuring the projected distance in 
Considering the large number of ionizing stars in LH 10 



Fig.|2 

(Sect. |2.2| ), it is surprising to find such a low ionization parame- 
ter. However, we note that the flux of all the ionized gas tracers 
(Ha, [O iii]88,opt. [N III], and [N ii]i22) is in fact only twice greater 
toward W2 than toward the surrounding regions, so half of the 
ionization must come from an extended gas component, which 
has likely a lower U value because of geometrical dilution. We 
thus expect a very low-density component that drives U lower 
(see also Sect.|7]). Choosing a density of rie = 25 cm"^ and a dis- 
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Fig. 8. The line flux ratios plotted against the ionization parameter U for radiation temperature of T^f^ = 40000 K (orange models) 
and 45000 K (blue models). The horizontal lines show the observed values toward specific regions illustrated in previous figures, 
while the associated vertical bars show the range of values within each region. The light gray horizontal rectangle shows the range 
of values observed across the entire N 1 IB nebula. The dark gray vertical rectangle indicates the range of U values agreeing the best 
with the observed line ratios (see text). The vertical arrow shows the eff'ect of dust extinction on the observations for Ay = 1. The 
Wl region was not covered by the [O iiiJopt HST/WFPC2 observation (Fig.|6]), and [N ii]i22 was only detected by the HerschelfPACS 
follow-up pointed observations (Fig. [4]). 



tance of 30 pc (size of the PACS map coverage), log U = -2 
corresponds to about 26 05 V sources, which is a realistic num- 
ber considering the O star census in N 1 IB (Sect.|7]). 

[N ii]i22/[N III] is another useful tracer of the physical con- 
ditions in the ionized gas, as it usually provides an estimate 
of Teff (e.g., Rubin et al. 1994). This ratio is also geometry- 
dependent because [N ii]i22 emits close to the ionization front. 



so that the nebula should be radiation-bounded for the derived 
physical conditions to be meaningful. We now test this hypothe- 
sis. [N ii]i22 was only detected toward the brightest regions Wl 
and W2 (Fig.|4]). Figure[8|) shows that [N ii]i22/[N iii] toward Wl 
and W2 is compatible with an ionization parameter in the range 
[-2.5,-1.8] (He ~ 100 cm"^ and T^^ = 40000 K). We obtain 
the same range from the [N iii]/[0 iii]88 diagnostic plot, imply- 
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ing that [N ii]i22/[N iii] provides a good estimate of the physical 
conditions toward Wl and W2; i.e., the nebula must be radiation- 
bounded toward these two positions. This agrees wi th th e opti- 
cal images that show multiple ionization fronts (Sect. 2J_). If the 
nebula was not radiation-bounded, we should have observed a 
lower [N ii]i22/[N iii] ratio, which in turn would have translated 
into a larger U parameter. Unfortunately, we cannot conclude 
anything about the other regions, since [N ii]i22 was detected 
only toward Wl and W2. 

The diagnostic plots in Fig. [8] give useful indications about 
the excitation conditions within the compact source Wl. The 
ratios [N ii]i22/[N iii] and [N iii]/[0 iiilgg toward Wl are simi- 
lar to the values within W2. Assuming that the same radiation 
field temperature holds for both regions (~ 45000 K), this im- 
plies roughly the same ionization parameter. Given the small size 
of Wl (1 pc; Fig.|2]), we could have expected a much larger U. 
That U is the same toward Wl and W2 is a strong indication 
that sources within Wl do not dominate the ionization of the 
region. Instead, it is likely that the ionization fronts are due to 
sources outside Wl. These sources are very likely the massive 
stars in the "footprint" cluster and in LH 10. This was hypothe- 
sized by Naze et al. (2001) based on the orientation of the ion- 
ization fronts as seen from the optical images. Sources within 
Wl must be much colder than 40000 K, resulting in a larger U 
parameter but a much weaker and, in fact, negligible [O mlgg 
flux. 



6.3. Evidence of obscured ionized gas 

We now investigate the dust extinction by comparing [O mlgg 
and [O III] opt- The FIR line is not aff'ected by dust extinction and 
potentially provides a better handle on the spatial distribution 
of the ionized gas. Furthermore, understanding the extinction 
across N IIB allows us to use the [O iii]opt/HQf ratio as a diag- 
nostic tool for the physical conditions. 

The relative variations of [O iiiJopt and [O miss, which intrin- 
sically depend on the electron density and temperature (Sect.[5|, 
is then modified by dust extinction, which significantly aff'ects 
the emission from the optical line emission. The variation in 
[O III] 88 and [O iiiJopt between LH 10 and W2 is remarkably dif- 
ferent, the FIR line emission being only 2.5 times lower toward 
W2 than toward LH 10, while the optical line is ten times weaker 
toward W2 than toward LH 10 (Figs. |4] and [6]). This could be due 
a priori to a variety of parameters (diff'erent line critical densities, 
diff'erent electron temperature, extinction by dust). 

Figure^ : shows that [O 

iii]88/opt is systematically underesti- 
mated by all models. The most relevant models for comparison 
to observations have densities around (~ 10 - 1 00 cm~^; long- 
and short-dashed lines) and U - -2 (Sect. |6.2| ). As shown in 
Fig. [7] [O iii]88/opt varies with the electron density and tempera- 
tureTTsamis et al. (2003) find a large discrepancy between opti- 
cal recombination vs. collisionally-excited lines of O^, but they 
note that this diff'erence cannot be explained either by tempera- 
ture fluctuations in a chemically and density homogenous nebula 
or by high-density clumps. 

Another caveat concerns the metal abundance. Although the 
observed [O iii]88/opt ratio could be reconciled with models by 
using a higher metallicity (since heavy elements provide more 
cooling, which in turn reduces the optical line emission), we 
consider this hypothesis as being unrealistic given the well con- 
strained abundances in N IIB (Tsamis et al. 2003). For illustra- 
tion, a solar abundance would be required for the models to reach 



Given the difl'erent wavelength ranges of [O iii]88 and 
[O III] opt, extinction by dust plays the most important role. The 
efl'ect of dust extinction with Ay = 1 is shown in Fig.jSj:. The 
average color excess towards stars across the NllB nebula is 
£'(B-V) ^ 0.2, while some stars have a significantly greater red- 
dening (Parker et al. 1992; Lee 1990). This can be compared to 
^(B-V) ^ 0.05 toward the central cluster LH 9 of N 1 1, which 
is essentially foreground galactic extinction (Parker et al. 1992). 
Tsamis et al. (2003) find that the B aimer ratios in the optical 
spectra (long slit of 5.6' x 1.5'') agree fairly well with the theoret- 
ical values in the case B with a Galactic foreground extinction. 
From the same Balmer ratios, we estimate that the extinction 
cannot be much more than E(B - V) < 0.15 in the area probed 
by the long slit. From this value of £'(B-V) = 0.15, we estimate 
that as much as 38% of the [O iiiJopt emission is extinguished 
by dust by using the extinction curve from Cardelli, Clayton 
& Mathis (1989) with a total- to- selective extinction Ry = 3.1 
(Howarth 1983). This leads to a lower [O iii]88/opt by a factor 
^ 1.6, which is high enough to reconcile most of the observa- 
tions with the models. 

Toward a more detailed analysis, we find that [O iii]88/opt 
ratio toward the stellar cluster LH 10 is a factor of 2 above 
any models. In particular, the observed ratios lie above the mod- 
els with the low densities expected in that region (Sect.[23). 
Assuming that the diff'erence between the models and the ob- 
servations is driven by interstellar dust extinction, such a factor 
translates into a visual extinction of Ay ^ 1 or E(B-Y) ^ 0.21. 
This determination is in good agreement with the average color 
excess £'(B-V) = 0.2 toward the stars. 

The situation toward W2 is more complex. If the line emis- 
sion is dominated by low-density ionized gas, then almost no 
dust extinction is required since the observed [O miss/opt is close 
to the 10 cm"^ track. On the other hand, if the emission is dom- 
inated by denser gas (Sect. 2.1 ), then a factor 1.5-2 could be 
required, translating into E(B-Y) ~ 0.12 - 0.2. Finally, the 
[O iii]88/opt ratio toward the northern region is also close to 
the lowest density tracks, requiring only little dust extinction 

(Ay ~ 0.5). 

We thus get the picture of a low-density nebula emission that 
is hardly attenuated by dust extinction, except toward the stellar 
cluster LH 10 and possibly W2. Is the required extinction com- 
patible with the amount of dust available? We now determine the 
maximum dust extinction across the region. The maximum ex- 
tinction is given by the total amount of dust available. We built 
the extinction map from the modeling of the spatial distribution 
of the dust mass derived by Hony et al. (in preparation). The 
grain emission is constrained by the Spitzer and Herschel data of 
Nil. The grain composition is the 'Amorphous carbon model" 
of Galliano et al. (2011), since it is the model that gives gas- 
to-dust mass ratios that are consistent with the elemental abun- 
dances in the LMC. In addition, it provides a conservative esti- 
mate, since it gives a lower dust mass, hence a lower Ay. The 
dust mass derived from SED modeling is turned into an Ay by 
using 



= ^ext(^) X ^dustMpix ^ 



(2) 



the observed [O iii] 



/opt 



toward Wl and LH 10. 



where A^^^ = 124 pc^ is the surface of the pixel used for the 

modeling, and A:ext(^) = 3175 m^/kg is the V band opacity of 
the grain mixture (Fig. A.l of Galliano et al. 201 1). Equation ^ 
corresponds to the maximum extinction since it corresponds to 
the Ay of the entire LMC toward the line of sight. The Ay^^ cal- 
culated from the dust modeling ranges from ~ 1 to ~ 16 in the 
area observed by the PACS FIR lines, with a median value of 3. 
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^max _ ^ 4 toward LH 10, ~ 5 - 16 toward W2, and -3-4 Table 5. O stars with known spectral types, 
toward the northern region. We call here that the actual extinc- 
tion Ay corresponds to the dust fraction between and the relevant 
emission source (here the emission-lines from the ionized gas). 
Without a precise knowledge of the morphology of the region, 
this fraction lies anywhere from to 1 . We nevertheless conclude 
that there is enough dust to explain an extinction of Ay ~ 1 to- 
ward LH 10 and W2. 

Correcting for dust extinction with Ay ~ 1 toward LH 10 
and W2, we find that [O iii]opt/Ha toward these two regions is 
in reality significantly greater than toward the northern region 
(Fig.[8]l). This is simply explained by the larger distance between 
the stellar cluster and the northern region, which reduces the ion- 
ization parameter. 



7. Spatial extent of the [O m] emission 

We observe that the [O iiilgg emission is spread throughout the 
Nil nebula. Is the [O iiilgg emission compatible with the number 
of known massive stars and their spatial distribution? The ISM 
density structure plays a dominant role in controlling the mean 
free path of FUV photons. The presence of [O iiilgg emission sig- 
nificantly far away from any massive star is a strong indication 
that photons are able to travel through a low-density material 
and/or that the ISM is significantly porous. We now attempt to 
model the [O iiilgg emission around each massive star in order to 
constrain an upper limit on the density above which FUV pho- 
tons are not able to travel to large distances. 

We consider the O stars identified by Parker et al. (1992) 
(Table|5]) since cooler stars are not able to ionize significant 
amounts of oxygen into O^^ (Sect.|2.2|). It must be noted 



Name"" 


Type^ 




f 


^([0 III]88) 






(pc) 




(xlO-^4 Wm-2) 


Eastern half^ 




22.6 (20.0) 


0.5 


24.4 


3173 (Nl 1-080) 


07 V 


7.2 (6.0) 


1.0 


1.2 


3224^ (Nl 1-050) 


04 V 


13.4(11.2) 


0.7 


4.7 


3204 (Nl 1-048) 


06.5 V 


11.0 (9.0) 


1. 


3.5 


3168 (Nl 1-032) 


0711 


12.2(11.3) 


1. 


6.2 


3223 (Nl 1-013) 


08 V 


5.8(4.1) 


1. 


0.2 


3209" 


03 III 


18.8 (16.2) 


0.7 


16.6 


Western half^ 




27.8 (24.5) 


0.4 


29.3 


3061 (Nl 1-031) 


0N21I1 


18.8(16.3) 


0.3 


6.5 


3058 (Nl 1-060) 


03 V 


17.4 (15.0) 


0.3 


6.6 


3053 (Nl 1-018) 


0611 


17.7 (15.0) 


0.3 


6.2 


3042 (Nl 1-087) 


09.5 V 


3.2 (2.2) 


0. 


0. 


3100 (Nl 1-038) 


05 111 


10.2 (8.2) 


0.5 


1.4 


3070 


06 V 


11.0 (9.0) 


0.3 


1.3 


3126 


06.5 V 


10.2 (8.2) 


0.7 


2.0 


3120" 


05.5 V 


12.0 (10.0) 


0.7 


3.8 


3073 


06.5 V 


10.2 (8.2) 


0.4 


1.2 


3089 


08 V 


6.8 (4.9) 


0.5 


0.2 


3102 


07 V 


9.1 (7.0) 


0.6 


1.0 


3103 


09.5 IV 


3.4 (2.3) 


0.9 


0.07 


3115 


09 V 


4.7 (3.3) 


0.9 


0.2 


3123 


08.5 V 


5.8(4.1) 


0.9 


0.2 



that unidentified embedded massive stars can contribute to the 
[O III] 88 emission. Furthermore, for several stars, the identifica- 
tion could be uncertain because of multiple stars combined in the 
same spectrum. While the spectral type derived from the com- 
posite spectrum can be correct, the number of sources (and thus 
the number of ionizing photons) might be diff'erent. The number 
of O stars given in Table[5] should be considered as a lower limit. 

We assume in the following that stars excite their own sur- 
rounding nebula (i.e., photoionization is radiation-bounded) and 
that the density is homogenous. The stellar properties (number 
of ionizing photons, temperature) are taken from Sternberg et al. 
(2003). Figure|9] shows the [O iii]88 emitting spheres around each 
O star in N 1 IB. The emissivity from the Cloudy models around 
each star was computed assuming a uniform spherical geometry 
and projected on the sky at the distance of N 1 IB. The final maps 
(Fig.|9| were created by coadding the individual models. 

It can be seen that models with a density of > 100 cm"^ pre- 
dict a small ionized volume as compared to the physical scale of 
N 1 IB. As the density decreases, FUV photons are able to travel 
farther away and ionize the ISM several tens of parsecs from the 
star. A density around 8 - 16 cm"^ seems to reproduce the phys- 
ical extent of N IIB well. The [O iii]88 extent for each O star 
is given in Tablejsjfor a density of 16 cm"^. The largest diame- 
ter at a density of 16 cm"^ reaches ^ 40 pc, for PGMW3061, 
3058, 3053, and 3209, which is close to the physical size of 
N 1 IB (^ 80 X 60 pc from the Ha image). We conclude that it is 
possible for the O stars to excite [O iii]88 throughout the PACS 
map, provided some low-density channels exist that are as low 
as < 16 cm"^. This value is significantly less than the densities 
from the optical doublet diagnostics (~ 100 cm"^. Sect. 2.3 ), and 



it suggests that the optical line measurements are biased toward 
the brightest and possibly densest regions. 



Spectral type and name are from Parker et al. (1992). The label in 
parentheses refers to Evans et al. (2006). When available, the spectral 
type was updated using Evans et al. (2006). Radius of the [O iiilgg 
emitting sphere containing 99% of the total [O iiilgg^ assuming an 
homogenous density of 16cm~^. The value between parentheses give 
the radius of the sphere containing 75% of the total [O iii]. Fraction 
of the [O III] 88 emission falling into the PACS map (assuming a density 
of 16cm~^). The model for each half of NllB considers all stars 
located in the same location, it is not a combination of the models from 
individual stars. Possible composite spectrum in Parker et al. (1992). 



Although the absorption of UV photons is mostly due to the 
gas rather than to the dust, dust does play an important role in 
determining the mean free path of UV photons. It is therefore in- 
teresting to consider a cluster of stars ionizing the ISM around it 
rather than a collection of separate individual stars. Models with 
stars gathered in clusters remove the assumption that each star 
excites their own surrounding nebula. We therefore now con- 
sider two main regions powered by their own stellar cluster. The 
western half of N 1 IB is dominated by the LH 10 southern clus- 
ter and the eastern half by the "footprint" cluster (Fig. [9]). The 
nebular density best agreeing with the physical size of N IIB 
is around < 16 cm"^. The radius corresponding to each cluster 
is reported in Table[5]and shown in Fig.|9] The volume ionized 
by the clusters is slightly smaller than the sum of the individ- 
ual volumes around each star for a given region (by 17% and 9% 
lower for the western and eastern halves, respectively). Although 
the exact location of the cluster centroid is quite uncertain since 
stars are in reality not co-spatial, the overall [O iii] extent of the 
two clusters encompasses the N IIB nebula relatively well. 

We conclude that massive stars in N 1 IB are able to produce 
[O III] 88 emission at large distances in the global low-density 
medium of NllB. Although higher density material does ex- 
ist, the low-density regions provide channels for the FUV pho- 
tons to permeate into the ISM. For instance, the northern region 
shows significant [O iii]88 emission, while there are no ionizing 
sources, implying that FUV photons were able to cross most of 
the N IIB nebula. The global density of < 16 cm"^ that we de- 
rive is an upper limit, since projection eff'ects prevent us from 
asserting the depth of the nebula in detail so that distances are 
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Fig. 9. [O iii]88 distribution modeled around each star with an homogeneous density ranging from 2cm"^ to 256 cm"^. The yellow 
polygon indicates the PACS coverage and the green polygon indicates the approximate extent of the optical body. The cyan circles 
in the n = 16 and 32 cm"^ panels represent the extent of the [O iii] -emitting zone when considering a cluster instead of the individual 
stars (see text). 



in fact lower limits to the actual distances. Slightly larger dis- 
tances would then require slightly lower densities in order to 
reach further. Another caveat should be added that we assumed 
a gas filling-factor of 1. If the gas were to fill only a fraction 
of the nebula, UV photons would travel even farther for a given 
density. 



8. Origin of the [C ii] emission 

Because of its low excitation temperature, [C ii] can potentially 
emit from several ISM phases (see introduction). The [C ii] emis- 
sion peaks toward W2, but it is dominated by extended emission 
(Sect.|5.2[). The [C ii] distribution is characterized by a sharp de- 



crease in flux toward the stellar cluster LH 10, contrasting with 
the flat [O III] 88 emission (Fig.|4]). In Fig.[TO| it can be seen that 
the spatial distributions of [O iii]88 and [C ii] are anti-correlated 
to first order, except toward W2 where all the FIR tracers are 
concentrated. The [O iii]88 emission thus appears to surround 
lower ionization regions, which might be associated to photodis- 
sociation regions. In the following we now investigate the possi- 
ble contribution from the ionized gas and from PDRs to the [C ii] 
emission. We use [N ii]i22 and [O iii]88 to trace the low-density 
ionized gas, and [O i]63 and PAH emission to trace PDRs. 
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Fig. 10. Three-color image of the following PACS lines, [C ii] 
(blue), [O i]63 (green), and [O miss (red). 



SA. Ionized gas component 

Most of the ionized gas in N 1 IB has a density around or lower 
than the critical density of [C ii] (80 cm"^ for collisions with e~; 
Table[T]), so we expect most of the in the ionized phase to 
de-excite radiatively in the ionized gas. Figure[TT] shows that the 
[C ii] emission is systematically underestimated by the models. 
We wish to call here that our models correspond to the photoion- 
ized gas in the H ii region an d tha t calculations were stopped at 
the ionization front (see Sect. 6^). An important warning is that 
[C ii] in the ionized gas is expected to emit close to the ion- 
ization front, so that our results strongly depend on the geom- 
etry and on the radiation-bounded nature of the nebula. While 
the compact regions Wl and W2 appear to be radiation-bounded 
(Sect. |6.2| ), the situation is uncertain for the more diffuse and ex- 
tended regions. Therefore, the model predictions for [C ii] emis- 
sion should be considered as upper limits except toward Wl and 
W2. With these caveats in mind, we assume a density between 
10 cm"^ and 100 cm"^, and from Fig.[lT]we estimate that the 
fraction of [C ii] originating in the ionized gas is ~ 15% toward 
LH 10 and ~ 4% toward the northern region, while it is < 4% 
toward Wl and < 2% toward W2. 

The [N ii]i22 line is a good tracer of the diffuse ionized 
gas from which [C ii] could also arise (e.g.; Oberst et al. 
2006). [N ii]i22 and [C ii] involve similar energies (ions exist- 
ing for energies > 14.5 eV and > 1 1.3 eV respectively), but they 
have significantly different critical densities in the ionized gas 

400 cm"^ and ^ 80 cm"^, respectively, for collisions with 
e~). Even for densities as low as ~ 50 cm"^, therefore, the de- 
excitation of [C ii] in the ionized gas is affected by collisions, 
and the [N ii]i22/[C ii] ratio thus increases with density until the 
critical density of [N ii]i22 is reached. [N ii]i22 was detected only 
toward the compact regions Wl and W2 with pointed obser- 
vations (Sect.[5]). The PACS footprint size unfortunately barely 
samples the extended emission around these two regions (Fig.|4]). 
FigurefTT] shows that [N ii]i22/[C ii] toward Wl and W2 is signif- 
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Fig. 11. The integrated emission of [C ii] (top) and the 
[N ii]i22/[C ii] ratio (bottom) are plotted against the ionization 
parameter U. Models show the expected values in the ionized 
gas of the H ii region. See Fig.[8]for the plot description. 



icantly lower than the models with Ue ~ 100 cm ^, confirming 
that most of [C ii] arises in PDRs toward these regions. The ISO 
observation (Sect.|T2]) shows that [N ii]i22/[C ii]< 0.03 in the 
80" beam, which includes W2. It thus seems that our results re- 
mains unchanged on a larger scale; i.e., [N ii]i22/[C ii] is lower 
than expected from the ionized gas, and most of [C ii] toward 
Wl and W2 originates in PDRs. 



16 



Lebouteiller et al.: Physical conditions in the gas phases of the H ii region LMC-N 1 IB unveiled by Herschel 



8.2. PDR component 

Our results suggest that most of the [C ii] emission in NllB 
does not originate in the ionized gas (Sect. |8.1| ). However, the 
uncertainty on the ionization structure and the low S/N ratio of 
the ionized gas tracer [N ii]i22 observation prevents a conclusive 
result for the most diffuse regions (i.e., other than the regions Wl 
and W2). We now investigate the relation between gas heating 
tracers (e.g., TIR and PAH emission) and gas cooling tracers 
(e.g., [C ii] and [O i]63) to understand the origin of [C ii] and 
also to constrain the dominant gas heating source in PDRs. 



8.2.1. Gas heating 

The gas heating in PDRs and in neutral atomic clouds is domi- 
nated by the photoelectric effect on dust grains (e.g., Hollenbach 
& Tielens 1999; Weingartner & Draine 2001), in which a pho- 
toelectron is ejected from the grain following the absorption of 
a FUV photon. The power absorbed by dust that goes into ISM 
heating, referred to as the photoelectric efficiency 6pe, is ~ 1%, 
and is the largest for small grains and neutrally charged grains 
(e.g., Tielens & Hollenbach 1985a). The gas heating is expected 
to originate about equally in PAHs (and PAH clusters) and in 
very small grains with radii < 100 A (Bakes & Tielens 1994). 

Although TIR emission mostly traces classical large grains 
in thermal equilibrium with the interstellar radiation field, the 
TIR emission is a good proxy for the gas heating rate in PDRs 
assuming a given photoelectric efficiency. The [C ii]/TIR ratio 
has often been used to approximate 6pe, under the assumptions 
that TIR traces the gas heating rate and that [C ii] dominates 
the gas cooling (e.g., Rubin et al. 2009 and Israel et al. 2012 in 
the LMC). The TIR map we used to estimate the dust extinc- 
tion in Sect. |6.3| is too coarse to compare to the FIR line maps 
(spatial resolution of 23'' as compared to ^ 11.5''; Fig.|4]). We 
therefore examined the correlation between TIR and single dust 
continuum bands in order to derive a higher resolution proxy 
for TIR. The following bands were considered: MIPS 24yum, 
PACS 100 yum and 160 yum, and SPIRE 250yum, with resolutions 
of 6'\ l'\ 12'', and 17" respectively. The continuum images 
were convolved to 23" resolution and resampled to match the 
TIR image. The correlation between TIR and each continuum 
band is shown in Fig. [12] We find that TIR correlates best with 
the lOOyum continuum on the scale of the entire Nil nebula 
and most importantly within the star-forming regions N 1 1 A and 
N 1 IB (Fig.[T2]). The scatter is also the smallest for TIR/100 yum. 
The linear regression fit shows that the flux density ratio is 
TIR/lOOyum = 2.70 ± 0.27. In the following, we use the 100 yum 
dust continuum image to trace TIR to a spatial resolution com- 
patible with the [C ii] and [O i]63 mapping observations shown 



in Fig. 12 



PAH molecules and PAH clusters are expected to be impor- 
tant contributors to gas heating. Rubin et al. (2009) find a good 
correlation between [C ii] and the Spitzer/IRAC 8.0 yum band 
throughout the LMC and proposes that PAHs dominate the pho- 
toelectric heating rate in PDRs. In the following, we revisit these 
findings by using the PAH map presented in Sect. 4.2 and by ex- 
amining both [C ii] and [O i]63 as potential coolants. 

8.2.2. Relation between dust heating and gas cooling 



Figure[T3] shows the correlation between gas cooling tracers 
([C ii], [O i]63, and their sum) and dust heating tracers (TIR, 
PAH). We use h ere the lOOyum dust continuum as a proxy for 
TIR (Sect.|8.2.T]). The ratio of cooling lines to dust heating trac- 
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Fig. 12. Top - TIR map of N 1 1 with a resolution of 23'' (see 
Sect. |6.3| ). Middle - Correlations between TIR and individual 
dust continuum bands. Black points correspond to positions in 
N 1 1 A and NllB and gray points to the rest of the nebula (see 
top panel). The dashed line shows the average value. Bottom- 
TIR map using the 100 yum dust continuum as a tracer. See Fig.|4] 
for the figure description. 
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ers provides a proxy for the fraction of the power absorbed by 
dust that is transferred into gas heating. The [O i]63/(TIR,PAH) 
ratio shows a large scatter (factor of 15 for [O iJes/TIR, ^ 5 
for [O i]63/PAH), with the lowest values toward the stellar cluster 



LH 10 and the northern region (Fig. 13 ), while the highest values 
are found in Wl and W2. The [C ii]/(TIR,PAH) ratio is signif- 
icantly tighter than [O i]63/(TIR,PAH), with a factor of ^ 6 for 
[C ii]/TIR and ^ 4 for [C ii]/PAH. Finally, the use of the sum 
[C ii]+[0 i]63 to trace the gas cooling provides the tightest rela- 
tion with both TIR and PAH emission, with a factor of 4 for 
([C ii]-F[0 i]63)/TIR and ^ 2 for ([C ii]-h[0 i]63)/PAH. We note 
that the relatively large scatter of ([C ii]+[0 i]63)/TIR is caused 



by the LH 10 region and its surroundings (see also Sect. |8.2.4| ). 
When ignoring LH 10, the scatter of ([C ii]+[0 i]63)/TIR is sim- 
ilar to that of ([C ii]-h[0 i]63)/PAH. 

The [C ii]/(TIR,PAH) and [O i]63/(TIR,PAH) ratios are 
anti-correlated in most regions. The ratios are also anti- 
correlated with respect to the TIR and PAH emission, with 
[C ii]/(TIR,PAH) decreasing with TIR and PAH emission, 
while the inverse is observed for [O i]63/(TIR,PAH). The anti- 
correlation between [C ii]/(TIR,PAH) and [O i]63/(TIR,PAH) is 
due to different PDR conditions, with the cooling being domi- 
nated either by [C ii] in the most diffuse regions or by [O i]63 in 
the densest regions and/or the regions with intense UV field (e.g., 
Kaufman et al. 2006). This result is compatible with the fact 
that [O i]63/(TIR,PAH) is relatively brighter toward Wl and W2, 
which are relatively dense regions located near massive stars 
(Sects. [2?T] |2.3| ). The cooling via [O i]63 represents ~ 25% of the 
total cooling rate across N 1 IB except toward Wl and W2 where 
it reaches ~ 50%. On large scales, we therefore expect [C ii] to 
be a fairly good proxy for the total cooling, which confirms the 
estimates made by Rubin et al. (2009). We keep in mind that 
the [O i]63 line could be affected by optical depth effects, more 
so than [O i]i45 and [C ii] (e.g., Abel et al. 2007 ). Ho wever, only 
low Ay values are found in N 1 IB (Ay < 1 ; Sect. |6.3| ), suggesting 
that optical depths effects are negligible at the spatial resolution 
probed by PACS (^ 9.5 - 13"). [C ii] and [O i]63 thus compen- 
sate for each other and, assuming these two lines are the domi- 
nant coolants at low Ay (e.g., Tielens et al. 1985a; Hollenbach 
& Tielens 1999), the sum [C ii]+[0 i]63 provides a proxy for the 
total gas cooling. 

8.2.3. PDR conditions 



The flat correlation between [C ii]+[0 i]63 and PAH (or, al- 
though to a lesser extent, TIR), together with the fact that Wl 
and W2 are dominated by PDRs (Sect.[8]T]), suggests that [C ii] 
and [O i]63 originate in PDRs throughout N IIB. The scatter of 
([C ii]+[0 i]63)/PAH is a factor of ^ 2, implying that the frac- 
tion of [C ii] originating in PDRs is at least ~ 50% throughout 
N IIB. The spatial anticorrelation between [C ii] and the warm 
photoionized gas probed by [O miss (Fig. 10 ) suggests that PDRs 
are mainly located within the ionized gas nebula (i.e., they are 
not located at the edge of the N 1 IB region). 

The small scatter of ([C ii]+[0 i]63)/PAH further implies that 
the photoelectric efficiency in PDRs is uniform within a factor 
of ^ 2. This implies that the grain charging parameter, related to 
the UV field intensity and to the electron density (e.g., Tielens 
& Hollenbach 1985a), is also uniform across the region. Rubin 
et al. (2009) find a bias toward lower [C ii]/TIR ratios in star- 
forming regions of the LMC, which they ascribes to relatively 
low photoelectric efficiency. The [C ii]/TIR in the entire Nil 
nebula found by the authors, ^ 0.56%, agrees with the highest 
values we measure in N 1 IB. It is therefore likely that the large 



beam used by Rubin et al. (2009), 14.9', was filled by the most 
diffuse regions where the cooling is dominated by [C ii]. Our 
results also imply that [C ii]/TIR alone (or [C ii]/PAH) does not 
trace the photoelectric efficiency well on small enough scales. 

We do not find any significant evidence of a lower photo- 
electric efficiency in the regions with the greatest expected UV 
field intensity, such as W2 or LH 10, as compared to more qui- 
escent regions, such as the northern region. We note, however, 
that the ionization parameter in most regions, including W2, 
indicates a dilution of the UV field on large physical scales 
(Sect. |6.2[ ). It is therefore possible that variations in the UV field 
intensity across N IIB are not large enough to modify signifi- 
cantly the grain charging parameter. From Fig. [13] the average 
ratios are ([C ii]-h[0 i]63)/PAH= 7% and ([C ii]-h[0 i]63)/TIR= 
0.55%, which can be considered as proxies for the photoelec- 
tric efficiency, if either PAHs or dust probed by TIR dominate 
the g as heat ing. We i nvestigate the dominant heating source in 
Sects.EallandESS 

8.2.4. LH10 

The case of the stellar cluster LH 10 is critical to understand- 
ing the dust component responsible for the gas heating in PDRs. 
LH 10 is dominated by warm photoionized gas, and the gas heat- 
ing is due to photoionization in addition to the photoelectric ef- 
fect. Most important, relatively warmer dust contributes to the 
TIR emission that is not associated with the heating of the gas 
in PDRs probed by [C ii] and [O i]63. In Sect. [0!2l we found 
that LH 10 drives the scatter on the ([C ii]+[0 i]63)/TIR ratio up. 
In contrast, the ([C ii]+[0 i]63)/PAH ratio shows a small scat- 
ter throughout N IIB (factor of ^ 2). This shows that the PAH 
emission describes the gas heating well even toward the stellar 
cluster where other dust components exist (in particular stochas- 
tically heated small grains and warm grains in thermal equilib- 
rium with the interstellar radiation field). 



The dust model presented in Sect. |6.3| indicates that the PAH 
abundance varies by a factor of ~ 10 across N IIB, with a min- 
imal value toward LH 10 where most PAHs are likely destroyed 
by the radiation field. Still, ([C ii]+[0 i]63)/PAH remains con- 
stant within a factor of 2 as compared to other regions in 
N IIB. PDRs therefore exist toward LH 10 and PAH emission 
is a good proxy for the gas-heating rate in these PDRs even to- 
ward such an extreme environment. The PAH/TIR ratio is the 
lowest toward LH 10, and only a fraction of TIR must originate 
in PDRs, abo ut - 4 0%. We investigate the dominant heating 
source in Sect. 18.2.5] 



8.2.5. Dominant gas heating source 

The gas hea ting sh ould be dominated by very small grains and 
PAHs (Sect. 8.2.1 ). The lAjim dust continuum, which is often 
used to trace warm dust grains and stochastically heated small 
grains, correlates poorly with [C ii]+[0 i]63, with a scatter ~ 5 
times larger than the correlation with PAH emission (not shown 
in Fig. 13). On the other hand, we find a remarkable correla- 
tion between [C ii]+[0 i]63 and PAH emission, suggesting that 
PAHs could dominate the gas heating over the dust components 
probed by TIR emission (Sect. 8.2.2 ). However, a similarly small 
scatter exists with TIR emission if we ignore the LH 10 region 
(Sect. |8.2.4| ). Croxall et al. (2012) also find a smaller scatter 
with ([C ii]-h[0 i]63)/PAH than with ([C ii]-h[0 i]63)/TIR and 
concludes that the integrated dust emission does not describe 
well the photoelectric heating where [C ii] and [O i]63 originate. 
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Fig. 13. Pixel-to-pixel correlation between the PACS cooling lines ([C ii], [O i]63, and their sum) and 2 quantities used to trace the 
gas heating, i.e., the TIR emission {left panel) and the PAH emission {right panel). The solid curve represents the linear regression 
for a given plot and the dashed curves the linear regression of the other ratios for comparison. The colors correspond to specific 
regions illustrated in Fig. [12] The symbol size is proportional to [C ii]+[0 i]63 {top), [C ii] {middle), and [O i]63 {bottom). The gray 
rectangles on the left side show the values observed by Croxall et al. (2012) in 2 super-solar metallicity nearby galaxies (see text). 



Instead, photoelectrons ejected from PAHs would dominate the 
gas heating (see also Helou et al. 2001; Rubin et al. 2009). We 
consider in the following PAHs and dust grains probed by TIR 
emission as potential gas-heating sources. 

The total gas-cooling rate Ag, traced by [C ii]+[0 i]63, is 
equal to the gas heating rate resulting from the photoelectric ef- 
fect on dust grains probed by TIR emission (F™) and on PAHs 

(rr): 



Ag-1 -hi - 6pE,TIR 1 d + ^PE,PAH 1 d , 



(3) 



where F^^^ and F^^^ are the dust heating rates for dust grains and 
PAHs, and 6pe,tir and 6pe,pah are the corresponding photoelec- 
tric efl&ciencies. Assuming that both dust components contribute 
to the gas heating, we conclude that upper limits on 6pe,tir and 
6PE,PAH are given by the ([C ii]+[0 i]63)/(TIR,PAH) ratio, i.e., 
0.55% and ^1% respectively. 
The PAH/TIR ratio varies by a factor of 5 across 
N 1 IB (Fig. [14]), enabling us to study the influence of PAH 
and TIR emission as gas-heating tracers. Figure 14 shows that 
([C ii]-h[0 i]63)/PAH and ([C ii]-h[0 i]63)/TIR are a smooth func- 
tion of the PAH/TIR ratio. The standard deviation around the 
([C ii]+[0 i]63)/PAH linear regression is remarkably small, with 
only 6%, while it is 15% for ([C ii]-h[0 i]63)/TIR. The top 
panel in Fig. 14 shows the trend of the ([C ii]+[0 i]63)/TIR and 
([C ii]-h[0 i]63)/PAH ratios versus PAH/TIR. When PAHs domi- 
nate the gas heating, the ([C ii]+[0 i]63)/TIR ratio increases with 
PAH/TIR, while the inverse is true for ([C ii]+[0 i]63)/PAH. We 



normalized each ratio to its average value and calculated that the 
linear regression slope is ^ -3.6 for ([C ii]+[0 i]63)/PAH, while 
it is ^ 10.3 for ([C ii]-h[0 i]63)/TIR. From this, we conclude that 
PAHs dominate the gas heating over the dust component probed 
by TIR emission. This implies that 6pe,pah is close to 7%, and 
therefore a factor of > 12 more than 6pe,tir. It is possible that 
^PE,TiR is larger if we would only a ccount for the gas heating 
from the smallest grains (Sect. 8.2. 1| ). 

Our results suggest that classical dust grains, which are dom- 
inating the TIR emission, do not dominate the gas heating. The 
fraction of the TIR emission corresponding to smaller grains 
certainly scales better with [C ii]+[0 i]63 than the total TIR. 
Unfortunately, in the lack of constraints on the abundance of 
small grains, we cannot conclude anything about the relative im- 
portance of PAHs and very small grains as gas-heating sources. 

9. Conclusions 

We present the HerschelfPACS spectroscopic observation of the 
H II region N 1 IB in the LMC. The spatial distribution of the FIR 
lines shows remarkable difl'erences: 

- The emission of [O miss, [N iii], and [C ii] is quite flat 
across the map, suggesting it is dominated by extended emis- 
sion. This contrasts with the [O i]63 emission, which is con- 
centrated in a few compact regions. [O i]63 also shows an 
extended component, about ten times fainter than the peak 
emission. 
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Fig. 14. Top - Expected trend between 

([C ii]+[0 i]63)/(TIR,PAH) against PAH/TIR, assuming that 
either PAHs or the dust grains traced by TIR dominate the gas 
heating. Middle and bottom - Pixel-to-pixel correlation be- 
tween ([C ii]-F[0 i]63)/(TIR,PAH) and PAH/TIR. The solid line 
represents the linear regression and the dotted line represents 
the median value. See Fig. [13] for the plot description. 



- The spatial distribution of [O nils 
correlated with that of [C ii]. 



appears to be anti- 



In the present paper, we investigate the extended emission 
of the FIR tracers, while in a second paper, we will investigate 
the compact regions and the PDR properties. The main results 
of the present study follow 



1. The [O III] 88 line, with its low critical density, is an excel- 
lent probe of the extended ionized gas. It is the brightest line 
across the region, in particular a factor of four brighter than 
[C ii] (in some regions 20 times brighter). 

2. By comparing [O iii]88 to the optical line [O iii] 5007 A, we 
find an extinction of Ay ~ 1 toward W2 and the stellar cluster 
LH 10, implying that a significant fraction of the ionized gas 
is hidden in the optical. The dust mass measured in N IIB 
is compatible with such extinctions. The extinction Ay is be- 
tween -0.5 and ~ 1 toward the other regions. 

3. The spatial extent of [O iii]88 emission could be reproduced 
by considering all the confirmed O stars, provided the den- 
sity of the ionized gas is generally lower than < 16 cm"^. 

4. We used the [N iii]/[0 iii]88 and [N ii]i22/[N iii] ratios to es- 
timate the physical conditions in the ionized gas. The aver- 
age ionization parameter is compatible with relatively low- 



density (< 25 cm"^) gas located on average < 30 pc away 
from the ionizing sources. 

5. By comparing [C ii] to the low-excitation, low-density ion- 
ized gas tracer [N ii]i22, we find that [C ii] mainly arises from 
PDRs across Nl IB. 

6. We examined the relation between the gas cooling in PDRs 
via [C ii] and [O i]63 and the gas heating traced by TIR emis- 
sion and PAH emission. We find a tight correlation between 
[C ii]+[0 i]63 and the PAH emission. Assuming [C ii] and 
[O i]63 are the dominant coolants at low Ay, the sum of [C ii] 
and [O i]63 provides a proxy for the total cooling, with [C ii] 
being the dominant coolant in most regions characterized by 
a diff'use ISM, and [O i]63 contributing significantly (up to 
^ 50%) in the brightest and densest regions. 

7. The scatter of the ([C ii]+[0 i]63)/TIR ratio is greater than 
that of ([C ii]+[0 i]63)/PAH, which is mostly due to the stel- 
lar cluster LH 10 where another dust component, unrelated 
to PDRs, drives TIR to higher values. 

8. We find that PAHs dominate the gas heating in PDRs, with 
a photoelectric efficiency of ~ 7%. 
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A.1. Transient correction 

Transients are cosmic ray hits to the detector that modify the 
pixel response over a given period of time (see the works of 
Coulais & Abergel 1999; 2000 for the ISOPHOT and ISOCAM 
detectors on ISO). Two types of transients are observed in the 
PACS data, "dippers" that result in a lower pixel response and 
"faders" that result in a larger pixel response. The response after 
a dipper hit normalizes after a short time (about the timescale of 
a line scan), while the response after a fader hit takes consider- 
ably longer to settle (near the observation itself or even more). A 
fader is often seen after the calibration block, which observes a 
bright source. The transient correction is only necessary for un- 
chopped scan observations since the fast chopping in chop/nod 
observations allows reliable removal of the transients (dippers 
and faders). 



A. 1.1. Faders 

Faders are responsible for large jumps in the flux, which can 
render a spectral pixel useless unless its signal is aligned with the 
others. HIPE includes a task called specLongTermTransient that 
minimizes the scatter of points by using a series of exponential 
functions to model the "faders" transients (Fadda & Jacobson 
2011). However, it is possible that several transients overlap in 
time, which makes it difficult to find the number of exponential 
functions required and the parameters of each function. This is 
the main reason we developed a special routine in P AC Smart that 
performs the transient correction on the signal timeline for each 
of the 16 spectral pixels and each of the 25 spaxels. 

P AC Smart includes a correction tool that uses a multiresolu- 
tion algorithm to straighten the signal (Starck et al. 1999). The 
smallest resolution scales are subtracted to remove the noise and 
the largest scale is subtracted to remove any fader-type signal 
(Fig. |A.l| ). The main advantage of multiresolution over an ex- 
ponential curve fitting is its ability to remove the signal from 
several faders that are overlapping in time. Furthermore, it pro- 
vides an unbiased approach since the faders do not need to be 
identified beforehand. 



A.I. 2. Dippers 

The presence of dippers dominates the S/N in the timeline (thus 
in the spectra as well). Dippers are searched for after each 
flagged pixel. An exponential curve is fitted with free parame- 
ters. If the timescale parameter is consistent with a dipper signal, 
the fit is used to correct the signal (Fig.|A.2). 



A.2. Line fitting 



Appendix A: The PACSman suite for Herschel/PACS 
spectroscopy data analysis 

PACSmar0 is an IDL package designed to provide an alterna- 
tive for several reduction and analysis steps performed in HIPE 
on PACS spectroscopic data. The following operations are cur- 
rently included: transient correction, line fitting, map projection, 
and map analysis. All of the observation modes are supported: 
wavelength switching (now decommissioned), unchopped scan, 
and chop/nod. 



^ Available at http://www.myravian.fr or by email contact (myra- 
vian @ gmail. com) . 



A3) is 



The input data for line fitting and map projection (Sect, i 

the Level 1 cube reduced with HIPE. At this level of reduction, 
the data is calibrated in flux and in wavelength. The bad pix- 
els are masked according to the HIPE reduction criteria (out- 
liers, glitches, saturated, noisy, etc.). The data cube contains one 
spectral window for each of the 25 spaxels and 16 spectral pix- 
els (and for a given raster position when mapping). The HIPE 
reduction corrects for the flat-fielding of the spectral pixels so 
that the combination of the spectral windows from the 16 spec- 
tral pixels provides one spectrum per spaxel. The data cloud for 
a given spaxel typically contains a few million points (for line 
spectroscopy). 

The line fitting in PACSman is performed on the full data 
cloud; i.e., the data is not rebinned beforehand. This is to keep 
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Signol vs, reodouts 




Fig.A.l. Example of fader transients occurring in unchopped 
scan observations. The signal is plotted against the readout num- 
ber for a given spaxel and spectral pixel. The irregular black 
curve is the unflagged data, while the flat green curve shows the 
corrected signal. Several transients are observed in this timeline, 
including the systematic transient after the calibration block on 
the extreme left. 



upper limits in faint lines. For a chop/nod observation, the two 
nods can either be combined in a single spectrum or considered 
individually. The fitting is performed as follows: (1) the contin- 
uum is fitted first around the line, and (2) the line -I- continuum 
profile is fitted using the continuum fit as a first guess. Example 
of fits are shown in Fig. [3] 

A. 3. Map projection 

For a mapping observation, the raster footprints (5x5 spaxel im- 
age of the line flux) have to be combined to create the final map. 
Maps are projected on a subpixel grid with a pixel size three 
times smaller (linear scale) than the spaxel size. The pixel size is 
about y\ This allows recovering the best spatial resolution pos- 
sible (see PACS ICC calibration document PICC-KL-TN-038). 

The grid orientation is chosen to be orthogonal with the 
raster maps so that the splitting of the spaxels requires as little 
rotation as possible. Since the spaxel positions and sizes are not 
regular, the orthogonality is only approximate. For a given pixel 
in the projected grid, the contributing fraction of each spaxel 
in the map is calculated assuming a uniform surface brightness 
within a spaxel (Fig. |A.3| ). Several determinations are thus ob- 
tained for each pixel in the projected grid. The final image sim- 
ply averages these determinations to obtain one value per pixel. 



Dip ^17, readout-8608, depth^7. 7474229 
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curve = 0.20721579, slope=0-0000O00, sdev = 0. JSS0e420, ch[= 1 .5758678 



A.4. Map analysis 

A tool is provided to analyze the spectral maps. Several tasks can 
be performed, such as aperture extractions with annulus subtrac- 
tion (see Figs. A.4 ). A 2D Gaussian surface can also be fitted to 
the data to estimate the source spatial FWHM. Finally, spatial 
profiles can also be plotted along a given cut. 



Fig. A.2. Example of a dipper fit. The profile is fitted by an ex- 
ponential curve whose parameters are compared to the expected 
behavior of a dipper (in particular the timescale). 



all the possible information on the data and on the masks during 
the line fitting process. The data cloud does not provide prop- 
agated errors for the individual measurements. We estimate the 
error by calculating the dispersion of measurements in a given 
wavelength bin (defined as 1/5 of the line FWHM). By default, 
a 5(T clipping is performed prior the line fitting. 

Line fitting is then performed with the MPFIT packag^ The 
continuum is adjusted first using spectral windows on each side 
of the line. The result is then used as an initial guess for the final 
fit that combines the continuum polynomial and the Gaussian 
curve. Default constraints exist that can be changed on the line 
FWHM, line central wavelength, and continuum range around 
the line. It is essential that constraints exist to provide reliable 

^ http://www.physics.wisc.edu/ craigm/idl/fitting.html 



22 



Lebouteiller et al.: Physical conditions in the gas phases of the H ii region LMC-N 1 IB unveiled by Herschel 




dx (■■) 



Fig. A.3. Example of a map projection. The observation (3x3 raster footprints) is shown with color symbols. The squares show 
the spaxel center (5x5 spaxels for one raster position), the large diamond shows the spaxels (1,1), and the dots show the center 
of the subspaxel (1/3 of the size of a spaxel). The information from different spaxels and rasters can be combined to populate the 
projected grid. Subpixels where there is no coverage are shown as gray squares. 




Fig. A.4. Example of a line cut across the map. The circles on the left image show the aperture used for extraction while the line 
shows the profile cut. 
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